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Some frustrated pyrochlore antiferromagnets, such as 
Y2M02O7, show a spin-freezing transition and magnetic ir- 
reversibilities below a temperature Tf similar to what is ob- 
served in randomly frustrated spin glasses. We present results 
of DC nonlinear magnetization measurements on Y2M02O7 
that provide strong evidence that there is an underlying ther- 
modynamic phase transition at Tf, which is characterized by 
critical exponents 7 ~ 2.8 and f3 ~ 0.8. These values are 
typical of those found in random spin glasses, despite the fact 
that the level of random disorder in Y2M02O7 is immeasur- 
ably small. 

*Present address: Department of Physics, University of 
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The past five years have seen a resurgence of significant 
interest devoted to the systematic study of geometrically 
frustrated antiferromagnets Geometric frustra- 

tion arises in materials containing antiferromagnetically- 
coupled magnetic moments which reside on geometrical 
units, such as triangles and tetrahedra, that inhibit the 
formation of a collinear magnetically-ordered state. The 
main motivation for the current interest in these systems 
stems from suggestions that (i) they may display criti- 
cal phenomena belonging to a "new" chiral universality 
class different from the universality classes of collinear 
magnets or (ii) the increased propensity of frus- 

trated antiferromagnets for quantum zero-temperature 
spin fluctuations compared to collinear antiferromag- 
nets might be sufficient to destroy Neel order and drive 
these systems into novel non-classical quantum disor- 
dered ground states 

Systems of classical Heisenberg spins residing on lat- 
tices of corner-sharing triangles or tetrahedra and anti- 
ferromagnetically coupled via nearest-neighbor exchange 
constitute particularly interesting cases of highly frus- 
trated antiferromagnets (see Fig. 1). Here, theory |3]|| 
and numerical work show that these systems do not 
order and remain in a "collective paramagnetic state" jj| 
down to zero temperature. Since, even for classical spins, 
these systems have such a small tendency to order, they 
are excellent candidates to display exotic quantum dis- 
ordered ground states ||,[|||. However, and perhaps 
most interestingly, experiments show that some nomi- 



nally perfect (i.e. disorder-free) [Q] pyrochlore antiferro- 
magnets @ exhibit a spin-freezing transition at some 
temperature Tf, below which they develop magnetic ir- 
reversibilities (see Fig. 1) and long-time magnetic relax- 
ation similar to what is found in conventional randomly 
frustrated spin glasses such as CuMn, EuSrS, and CdM- 
nTe 



111 



.0155 



.0145 



\ 
3 



.0125 



.0115 - 



II 



.0105 




5 10 15 20 25 30 35 40 

T(K) 




FIG. 1. The top figure shows the field-cooled (FC, squares) 
and zero field-cooled (ZFC, triangles) magnetization for 
Y2M02O7 in a magnetic field lOOOe. The bottom figure shows 
the three-dimensional network of corner-sharing tetrahedra 
formed by the Mo 4+ sublattice in Y2M02O7. 



Two important questions arise: Firstly, what is the 
microscopic origin of the glassy behavior in pyrochlore 
antiferromagnets? Is it due to the yet undetected mi- 
croscopic disorder inherent to any real material, or is it 
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intrinsic to the idealized perfect material? Secondly, ir- 
respective of the origin of the glassy behavior, one would 
like to know if the spin-freezing is strictly dynamical (i.e. 
where the system's relaxation time exceeds the time scale 
set by the experimental probe), or is it due to an under- 
lying thermodynamic transition characterized by a truly 
divergent (spin-glass) correlation length and time scale 
as is believed to occur in conventional disordered spin 
glasses pH? To address these two questions, we have 
measured the nonlinear magnetic susceptibility, Xni, of 
the pyrochlore antiferromagnet Y2M02O7. Briefly, we 
have found (i) that the freezing at Tf is well character- 
ized as a thermodynamic transition displaying a power- 
law divergence of the nonlinear susceptibility coefficient 
% 3 (T) ~ (T — Tf)-' 1 with 7 ps 2.8, and that the net non- 
linear susceptibility, Xni, exhibits critical behavior and 
temperature-field scaling properties close to Tf which 
gives a critical exponent, 8 ps 0.8. The values we ob- 
tain for 7 and 8 are typical of those found in conven- 
tional disordered spin glasses, despite any obvious mi- 
croscopic disorder in Y2M02O7 M. High temperature 
measurements of the susceptibility in the Curie- Weiss 

regime, X = C/{T - 9) show that 6 200K, with 

a Curie constant C giving an effective magnetic mo- 
ment of the order of 2.5 Bohr magneton per Mo 4+ Jl2| . 
This is very close to the theoretically expected value of 
2.8 Bohr magneton per Mo 4+ . We have recently per- 
formed high temperature susceptibility measurements in 
the range 10K— 800K which are in agreement with Hu- 
bert's results Clearly, these indicate that the sus- 
ceptibility, and thereby the spin-glass-behavior observed 
in Y2M02O7 is not due to a small fraction of the spins, 
but arises from the very large majority of the moments. 
Essentially, \ is small close to Tf — 23K because 9 is so 
large and since Y2M02O7 is an antiferromagnet (i.e 9 is 
negative). 

Y2M02O7 is a narrow band gap semiconductor where 
the Mo 4+ ions are magnetic, with an antiferromagnetic 
nearest-neighbor Mo-Mo superexchange, while Y 3+ is 
diamagnetic. Figure 1 shows the cubic unit cell of the 
Mo 4+ magnetic lattice where there is a magnetic Mo 
moment on the vertices of the tetrahedra. The 270mg 
powder sample of Y2M02O7 was prepared as described 
in Reference B. Neutron and X-ray powder diffrac- 
tion studies show that there is no measurable amount 
of oxygen vacancies or intermixing between the Y 3+ and 
the Mo 4+ sublattices ||. The random disorder in that 
material, most likely oxygen vacancies, is therefore be- 
low the 1% detectability level. The magnetization was 
measured using a commercial SQUID (Quantum Design, 
San Diego) magnetometer. The bulk magnetization of 
Y2M02O7 becomes history dependent below Tf ps 22K: 
the field-cooled (FC) and zero field-cooled (ZFC) mag- 
netizations measured in a field of lOOOe show a sharp 
breakaway below 22K (see Fig. 1) ||. 

To determine whether or not a true thermodynamic 
spin-freezing transition occurs around Tf ps 22K in 
Y2M02O7, we have measured the nonlinear susceptibility 



coefficient, \3{T), which is expected to show a power-law 
critical divergence at Tf fTl)] : 

X3 oc r" 7 , (I) 

with r = T/Tf — 1 and 7 > 0. \3 is extracted from the 
temperature, T, and field, H, dependence of the magne- 
tization, M(T,H) = X i(T)H - X3 (T)H 3 + X5 (T)H 5 - 
where Xi(^) is the linear susceptibility. Hence, the tem- 
perature dependence of X3(T) allows a determination of 
Tf and 7 In fact, all the nonlinear terms X2n+i with 
n > 1 must diverge at Tf, since both M(T, H) and H are 
finite quantities. It is therefore convenient to define a net 
nonlinear susceptibility, Xni, as 

it m - 1 M (T,H) 
Xn\{T,H) = 1 — . (2) 

Right at Tf , X ni has a power law dependence on H : 

Xm(T = T fl H)~H 2 / s , (3) 

where S is a second independent static critical exponent 
characterizing the spin- freezing transition |Tl| ]. Finally, 
and perhaps the most relevant test ascribing critical be- 
havior to the spin-freezing transition is obtained by seek- 
ing a scaling behavior of Xni of the form 

Xnl (T,H)o,H^ 5 T(r^^ 2 /H) , (4) 

where 8 is the spin-glass order parameter critical expo- 
nent p]J| . Here, J-(x) is the scaling function which must 
obey the following asymptotic behavior 

J-{x) = constant, x — > , (5) 
T{x) oca;" 27/(7+/3) , x -> 00 (6) 

in order that the scaling behavior has "physical con- 
tent" @, and that Eqs. | and | are recovered 0,0, 
hence giving the scaling relation 

* = 1 + 7/0 . (7) 

The magnetization data were obtained under field- 
cooling conditions. The field H, in the range [100 — 
7000] Oe, was switched on at high temperature (70K 
~ 37/ ), and kept constant during subsequent slow cool- 
ing at a rate 5 mK/s in a fixed field, and down to the 
temperature of interest. It took over an hour to go from 
70K to any of the (H — T) data point presented in our 
paper. This means almost twelve hours to get to the 
lowest temperature. The reason we cooled down in fixed 
field instead of working along conventional isotherms is 
to eliminate any possibility of a magnetic field hystere- 
sis of the superconducting magnet in the magnetometer 
affecting the results; the magnetic field was never de- 
creased during the whole experiment. Because of the 
irreversible and time-dependent nature of the system's 
response below Tf , only results in the temperature range 
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Tf < T < 37/ are included. Our results on the low 
temperature dynamical relaxation of the magnetization 
in Y2M02O7 will be reported elsewhere. Three consecu- 
tive cooling runs for fixed field were performed, with the 
magnetization data averaged over the three runs. 

Prior to doing any analysis, we have to deal with 
the fact that the interactions do not perfectly average 
to zero, as evidenced by Xi(^) n °t having a simple 
Xi ~ 1/T Curie law. It has been argued ]lT| , |l5[ sev- 
eral times that the leading (first-order) corrections to 
scaling coming from this non-zero averaging of the in- 
teractions can be eliminated by fitting X (T, H) to pow- 
ers of a 2n +i(T)xi[xiH} 2n for n > 0, instead of sim- 
ply X2n+iH 2n , and considering the critical behavior of 
o-2n+i{T) instead of X2n+i(T). For each temperature 
the field dependence of x &t small field was fitted with 
X = Xi - X3H 2 , giving a 3 = Xa/xl, and varying the 
upper limit of the field range to determine the limit of 
validity of this restricted fit beyond which higher X2n+i 
(n > 1) corrections become significant. The quality of 
our magnetization data did not allow us to determine 
the 05 = Xs/Xi coefficient with precision better than 
[50 — 100]%, and thus 05 data are not included here. 
We show in Fig. (2) the net nonlinear susceptibility, 
Xni(T,H), as defined in Eq. ||, with xi extracted from 
the fit M(T,H)/H = X i ~ X3H 2 , as a function of H 2 
for few temperatures above Tf. These results emphasize 
the large increase of the nonlinear susceptibility upon ap- 
proching Tf . We also notice that Xni is only linear in H 2 
up to a maximum field H+(T) whose value is rapidly 
moving to zero upon approching Tf due to the turning 
on of the X2n+i (n > 1) corrections which themselves 
diverge at T f as r (/3-«[7+/3]) f[J|. 
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FIG. 2. Net nonlinear susceptibility, Xnl = 1 — M/ X iH vs 
H 2 for the six temperatures indicated. 



Figure 3 shows a log-log plot of X3/X1 vs T/Tf — 1 for 
five different choices of Tf. The steepness of the curves 
increases as the chosen value for Tf is decreased. Choos- 
ing Tf > 23.0 would give us a curve with essentially no 
extended range of power-law behavior (X3/X1 x T ~ 7 )- 
As can be seen in the top of Fig. 1, values of Tf > 23. OK 
are definitely above the FC-ZFC break-away point where 
the magnetization data acquisition runs are reversible at 
the lowest field. Values of Tf < 21.0 are also not pos- 
sible to justify as they are clearly well below the FC- 
ZFC break-away point. Also, one observes an obvious 
upward curvature at the far left of the data with the 
choice Tf — 21K. Consequently, the log-log plot of the 
X3/X1 data strongly suggests that Tf is finite and be- 
tween 21. OK and 23. OK, and the extracted value of 7 
depends on the choice of Tf. To better quantify this, 
we have fitted X3/X1 00 (1 ~ T/Tf) 1 for a fixed range 
of value X3/X1 £ [2 — 30] for a fixed number of data 
points (solid line fits on Fig. 3). The goodness-of-fit, 
X 2 = Ei[« c 3 alc (n) - ar as -(n)} 2 , with a 3 = xz/xl and 
T} = 1 — Ti/Tf vs the choice of Tf is shown in the inset 
of Fig. 3, along with the corresponding variation of the 
critical exponent 7. As can be seen, a pronounced mim- 
imum in X 2 is seen as a function of Tf, and the best fit 
is obtained for Tf rj 22.2K , which gives a value 7 w 2.8. 
Overall, a conservative estimate gives Tf € [21.7— 22.7]K, 
giving 7 e [2.4 - 3.4]. 
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FIG. 3. Log-log plot showing the temperature dependence 
of az = X'ilX\ vs r = T/Tf — 1 for five different choices of 
Tf. T f = 21.0K, 21. 5K, 22. OK, 22. 5K, and 23. OK. The inset 
shows the goodness of fit parameters, X 2 , (filled squares) for 
the solid line fits in the main panel, and the exponent 7 (open 
circles) vs the chosen value for Tf (see comments in text). 



The power-law divergence of 03 
t < 0.05 for a choice of T f G [22.0 - 



= xa/xi saturates for 
23.0]K. A reason for 
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this levelling off of 0,3 is that the range of dominance of 
the term X3(T)H 2 to Xni falls below the smallest field, 
Hmin (-Hmin ~ lOOOe), for which good quality data were 
obtained. The increasingly important diverging higher 
order terms of alternating signs (\5, X7> etc) contributing 
to Xni then cause 03 to be underestimated when H + (T) 
becomes less than or equal to rs i? m in- Also, the slow but 
finite cooling rate inhibits the correct equilibrium value 
of Xni from being attained and this effect may be com- 
pounded with the previous one to produce a saturation of 
a 3 for t < 0.05. Overall, the observed behavior of (13 for 
Y2M02O7 seen in Fig. 3, as well as the uncertainty on the 
value of 7 are typical of what is observed in conventional, 
chemically disordered spin glasses. 

We now attempt to verify that the spin-freezing in 
Y2M02O7 is a legitimate critical phenomenon by seek- 
ing a data collapse and scaling behavior of the net non- 
linear susceptibility of the form given in Eq. § Choos- 
ing Tf = 22. 2K and 7 = 2.8, as found in the inset of 
Fig. 3, we can find a reasonable data collapse (scal- 
ing) of Xni with a choice of (3 = 0.75 ±0.10. Had we 
worked with other choices of T f G [21.7 - 22.7]K, and 
7 G [2.4 - 3.4], we would have found (3 G [0.6 - 0.9]. As 
found for 7 above, such a value for (3 is typical of that 
found in conventional disordered spin glases |11|. The 



scatter of the data at large x, 



-(7+/3)/2 



/H, arises be- 



cause these data are those that correspond to relatively 
high temperature (well above Tf) and small fields where 
the nonlinear magnetization is small and the experimen- 
tal error is the largest. In the limit x — > 0, we observe 
that T(x) approaches a constant value, hence confirming 
that Xni oc H 2 / 5 at T — Tf with Eq. @ obeyed. Tak- 
ing 7 = 2.8 and (3 = 0.75 at T = T f = 22. 2K, Eq. g 
predicts a value S ss 4.73. The inset of Fig. 4 shows a 
log-log plot of Xni{T = T f = 22.2)K vs H with this value 
of S = 4.73 (solid line). The fit is very good, confirm- 
ing that the collapse of the Xn\(T, H) data has physical 
meaning underlying a critical phenomenon which obeys 
Eq. 0. Also, in the limit of large x (i.e. at small fields), we 
observe that the asymptotic behavior of T{x) is consis- 
tent with the power-law T{x) oc x~ 2 7/(7+/3) (dashed-line 
in main panel of Fig. 4) as, again, it should be for the 
data collapse to have physical meaning and such that 
t-~<H 2 for T > T f and H -> 0. 
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FIG. 4. Nonlinear magnetization analyzed according to a 
scaling model for a nonzero spin-glass transition temperature 
with choices Tf =22. 2K, 7 = 2.8 and f3 = 0.75. The inset 
shows the log-log plot of xm(T = T f ,H) oc H 2/s with the 
value of S = 4.73 predicted from the scaling relation ^, with 
the values of 7 = 2.8 and (3 = 0.75. 



It is interesting to compare our results for the 
Y2M02O7 pyrochlore with those for the SrCrgGa^ig 
kagome system (SCGO) |l4|,[l(|, and for the site- 
ordered gadolinium gallium garnet magnet GdsGasO^ 
(GGG) @. Ramirez et al. |§ found a power-law di- 
vergence of X3 m SCGO with 7 w 2.4, while Martinez et 
al. |1J] recently argued that the freezing at Tf w 3.5K in 
SCGO is not associated with a divergence of X3_iX3 was 
found to increase by a factor 5 or so in Ref. (1J), and 
that this material does not exhibit conventional spin glass 
behavior. Also, contrary to what we find for Y2M02O7, 
Martinez et al. argued that the data-collapse (i.e. scal- 
ing behavior) they obtained for the net nonlinear suscep- 
tibility of SCGO was "unphysical" since the asymptotic 
behavior of their scaling function was inconsistent with 
what it should have been according to the scaling relation 
S = I + 7//3 Schiffer et al. (l7) found a large increase 
of X3 m GGG (6 orders of magnitude between 0.2K and 
5K), which they ascribe to a spin glass transition. How- 
ever, the temperature dependence of X3 in GGG is quali- 
tatively different than what is found in conventional spin 
glasses since X3 nas two maxima in GGG, while it is 
a monotonic function of the temperature in conventional 
spin glasses as well as here in Y2M02O7. Hence, from the 
point of view of nonlinear susceptibility measurements, it 
therefore appears that the spin glass behavior observed 
in Y2M02O7 resembles more closely what is found in con- 
ventional spin glasses than that which has been found in 
other topologically frustrated antiferromagnets, such as 
SCGO and GGG. 
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In conclusion, we have measured the nonlinear DC sus- 
ceptibility of the pyrochlore antiferromagnet Y2M02O7 
close to and above the freezing temperature, Tf ~ 22K, 
where this material shows spin-glass behavior. Our 
results show that the freezing transition observed in 
Y2M02O7 is well characterized by a power law diver- 
gence of the X3(T) nonlinear susceptibility coefficient, 
Xa(T) oc (1 - T/T/)-t, with a value 7 ~ 3.0 ± 0.5. 
This implies an underlying thermodynamic glass phase 
transition around 22K in this material. This is fur- 
ther supported by recent muon spin relaxation measure- 
ments where a large increase of the 1/T± muon depolar- 
ization rate at T as 20K has been observed, and which 
indicate a dramatic critical slowing down of the Mo 4+ 
moments fis| . The net nonlinear susceptibility data, 
Xni{T > Tf,H), can be collapsed onto a scaling func- 
tion from which we can extract the order pameter critical 
exponent (3 ~ 0.8 ± 0.2. Right at Tf, we find a behav- 
ior Xni ~ H 2 / 5 with a value of 8 ~ 4.7, which satisfies 
the scaling relation |?]. The values we find for 7 and f3 
are typical of what is found in conventional chemically 
disordered spin glasses jO]. 

We hope that our results will stimulate further experi- 
mental and theoretical studies to clarify the origin of the 
spin-freezing transition and the low-temperature glassi- 
ness seen in Y2M02O7 and other geometrically frustrated 
magnetic systems. 
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